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Abstract 

The optical vortex coronagraph is potentially a remarkably effec- 
tive device, at least for an ideal unobstructed telescope. Most ground- 
based telescopes however suffer from central obscuration and also have 
to operate through the aberrations of the turbulent atmosphere. This 
note analyzes the performance of the optical vortex in these circum- 
stances and compares to some other designs, showing that it performs 
similarly in this situation. There is a large class of coronagraphs of 
this general type, and choosing between them in particular applica- 
tions depends on details of performance at small off-axis distances and 
uniformity of response in the focal plane. Issues of manufacturability 
to the necessary tolerances are also likely to be important. 



1 Introduction 

Many areas of astrophysics need a coronagraph of very high dynamic range, 
to detect faint sources near bright ones. Examples include exoplanets, pro- 
toplanetary disks, and the structure of active galactic nuclei. These appli- 
cations cover a wide range of acceptable dynamic range, with planets being 
the most challenging. 
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A wide variety of coronagraphs has been proposed |T] and one particularly 
impressive one is the optical vortex coronagraph (OVC). This a remarkable 
example of diffraction at work [2[ |3l H] . A "vortex phase mask" ( VPM) is 
introduced at the focus at the focus of a circularly symmetric optical system. 
In the best-studied examples, this mask introduces a phase proportional to 
md, where 6 is an azimuthal co-ordinate in the focal plane, and m is an even 
integer called the topological charge. If the diffraction pattern at the focus is 
an Airy function, then one finds exactly, i.e. analytically, that all of the light 
in this diffraction pattern is diffracted outside the re-imaged pupil. A simple 
Lyot stop will thus extinguish all of the light from an on-axis source, at least 
for perfect optics. Actual implementations of a VPM will have imperfections, 
including possibly a loss of amplitude close to the vortex core, depending on 
how the VPM is implemented. Here I concentrate on an ideal system. 

The vortex phase mask appears to form the basis of a very elegant coro- 
nagraph. It has high performance by comparison with other designs [5j , and 
can be designed to be achromatic [6J and insensitive to aberrations. [7J Its 
conception can be seen as a continuation of related designs - the phase-mask 
coronagraph [8J and the four-quadrant phase-mask coronagraph. [9J 

It is of interest to examine the performance of the OVC on ground-based 
telescopes. These generally have two features which are not favourable for the 
OVC - a central obstruction, and images blurred by atmospheric turbulence. 
I will examine both aspects and show that the OVC behaves similarly to 
some other designs, but with good performance and some advantages even 
in these non-ideal circumstances. 

The non-ideal PSF — the central obstruction 

For any coronagraph, the loss of Strehl ratio from a central obstruction is a 
disadvantage. The OVC is no exception; in the reimaged pupil, the secondary 
is surrounded by a halo of light, whose total power is comparable to the power 
obstructed by the secondary. This result is illustrated by direct computation 
in Figure m 

This result can be examined analytically, as is discussed elsewhere in 
detail by Mawet et a/. [10]. Only the central obstruction need be added to 
their comprehensive treatment. I discuss a mask of charge 2 as an example. 
Using standard Fraunhofer diffraction, the amplitude, in a reimaged pupil 
plane beyond the VPM, is given by a generalized Hankel transform: 
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Figure 1: The intensity pattern in the reimaged pupil plane beyond a VPM, 
for an obstructed and unobstructed pupil. This was computed by direct 
FFT for perfect optics and charge 2 (top left), for charge 2 with a central 
obstruction (top right), for a quadrant phase mask (bottom left) and an 
optimum prolate apodization (bottom right). 



Air') = Ke^'^ / Tik) J2(kr')k dk 
Jo 

in which r and k are suitably scaled radial co-ordinates in the reimaged pupil 
and the preceding focal plane, and (3 is an azimuthal angle in the reimaged 
pupil. Constants irrelevant to the discussion are absorbed into K. T is the 
amplitude in the focal plane for an on-axis source. If this is an Airy function 
for an input pupil of radius i?, the integral becomes 

/•oo 

Air) = Ke^'^R / JAkR) Ukr') dk. 
Jo 
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(The second-order Bessel function arises because the amphtude introduced by 
the VPM has been taken to be e^*^ in this example. This topological charge 
of 2, via the angular integral of the Fourier transform, translates directly into 
the order of the Bessel function. The needed result is 





(see equation [5]). 

The infinite integral is a case of the Weber-Schafheitlin integral [11], and 

so 

A{r') = < r' < i? 

(1) 

= Ke^'f^i^f r'>R 

for the amplitude in the reimaged pupil. This is the result that motivates 
the optical vortex coronagraph; similar results can be obtained for other even 
topological charges of the VPM. The power diffracted outside R is \K\'^ttR'^, 
which is exactly equal to the power in the focal plane diffraction pattern 
{Ji{kR)/kRy. In the case of a centrally obstructed pupil, however, an on- 
axis source at infinity produces an amplitude in the focal plane that is given 
by 

^^^^ _ ^^RJi{kR) aJi{ka)^ 



k k ' 

in which a is the radius of the central obstruction. It follows that the ampli- 
tude in the reimaged pupil must be 

A{r') = < r' < a 

= -Ke^'^{^f a>r'<R (2) 

= Ke^^P U^f - {f,)A R>r' 



The total diffracted power is now \K\'^tt{R'^ — a'^), and the power diffracted 
into the reimaged pupil between a and R is |i^'p7ra^(l — (a/i?)^). The total 
energy diffracted into the halo of the secondary is nearly equal to the energy 
incident upon it, particularly for typical telescopes where a is much less than 
R. 
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A first-order approach to evaluating this issue is simply to examine the 
total power which passes through a Lyot stop. In fact, one might also consider 
blocking some of the diffracted light with an enlarged stop at the secondary; 
while this costs light, it might yield a net benefit in terms of attainable 
dynamic range. 



-1.5 

-2 

>< -2.5 
o 

-3 
-3.5 




0.2 0.4 0.6 

Light fraction used 



0.8 



Figure 2: The normalized diffracted power Xi plotted against the fraction 
of the light that is admitted through a centrally-obstructed Lyot stop. The 
curves are for linear central obstructions a/R of 0.2 and 0.1. The full line is 
for a VPM of charge 2, the dashed line for charge 4. 



Using Equation [21 the fraction of the power from an on-axis source that 
appears between p > a and R is 



X 




(3) 



This is the fraction of the energy of the the on-axis source that passes the 
proposed Lyot stop. This may usefully be compared to the fraction of the 
available light that is used in this scheme, which is 
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as is shown in Figure [21 For modest loss of light (~ 20%) the attainable 
dynamic range is about 100; this may be useful, especially if the dynamic 
range is more fundamentally limited by adaptive optics considerations. Los- 
ing light is not always a disaster if another parameter is improved, in this 
case dynamic range, but it is clear that the sacrifice would have to be very 
substantial to attain high dynamic range. 
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Figure 3: Calculations of the final point-spread function, for an on-axis ob- 
ject. The intensities are normalized to the diffraction-limited case with no 
central obstruction. Colours code the no- vortex case m = (black), charge 
2 (red), and m = 6 (green), all with a 20% linear central obstruction. 

A fuller evaluation looks at the point-spread function of an on-axis source, 
imaged through the VPM and a standard Lyot stop. This involves a Hankel 
transform of the amplitude in a > r' < R, as given by Equation [2] (or its 
extension for higher charges than 2). This integral is again of the Weber- 
Schafheitlin form. For charge 2, for example, the final focal-plane amplitude, 
in terms of a radial co-ordinate k', is proportional to 




Jijk'a) _ Ji{k'R) \ 
k'a k'R ) ' 



Examples of the PSFs are shown in Figure [3l where it is seen that the 
effects of the central obstruction are felt in different parts of the image, 
depending on the charge on the VPM. Depending on application, this level 
of light pollution may be acceptable. 
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For a ground-based telescope, partial correction for atmospheric turbu- 
lence may be a more relevant limitation to consider. Anticipating the dis- 
cussion of Section [2l Figure H] shows that the effect of a central obstruction 
could be quite perceptible. In this example the available "dark" region is 
completely filled by diffracted light. 





Figure 4: Simulation results, described later, for the on-axis long-exposure 
point spread function in a well-corrected AO case. The intensities are normal- 
ized to the diffraction-limited case. Colours code the no-vortex case m = 
(black), charge 6 VPM (solid green), and a quadrant phase mask (solid or- 
ange). Dashed lines are for with a 20% linear central obstruction. 
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Several methods are possible to recover the performance of the vortex 
coronagraph. For example, one might use a focal plane occulting stop as a 
second stage. Because the images formed on-axis through a VPM are quite 
wide (Figure H]) this stop would have to be large. 

Another option is to eliminate the central obstruction of the telescope, 
perhaps by using an off-axis Herschel-like system [T^. as is proposed for the 
Terrestrial Planet Finder mission. Large off-axis paraboloids are now be- 
ing manufactured for the Giant Magellan Telescope [13], which will contain 
six adjacent mirrors in an outer ring. From the Fourier shift theorem, the 
amplitudes in the reimaged pupils will simply be displaced versions of Equa- 
tion [H so that the issue will become the light pollution between mirrors 
(Figure E]). By using some subset of the available mirrors in a dilute aper- 
ture configuration, light pollution may be substantially reduced, although 
at the cost of losing light-gathering power. This is an example of how the 
central-obstruction problem may be handled with a more dilute aperture. 
Detailed calculation would be necessary to establish that the light pollution 
was acceptable for a particular configuration of apertures. An additional op- 
tion would be to include stops close to the focal plane that prevented light 
pollution from one aperture to the other, although this would limit field of 
view. 

The most attractive approach seems to be Guyon's [H] technique of 
phase-induced amplitude apodization. This reshapes the illumination of 
the pupil, using the complex aspheric surfaces which it is now possible to 
manufacture. Guyon gives examples in which the central obstruction os a 
two-mirror telescope is eliminated completely. Guyon was also concerned to 
tailor the resulting point-spread function so that it had low sidelobes, but 
there seems no reason why his technique should not be used in a simpler way 
to produce a pure Airy disk for an on-axis source. Apart from manufactura- 
bility, which is still costly, the main limitation to the method is the narrow 
field of view. If the pupil were reshaped to eliminate the central obstruc- 
tion, then an off-axis source will have a phase discontinuity at the centre of 
the reshaped pupil. This will limit the field to an angular extent given by 
the diffraction limit of a pupil of equal size to the obstruction. This will 
probably not be a limitation however in the very narrow-field applications 
of high-contrast coronography. It is also possible to avoid this field-of-view 
limitation by a second stage of pupil reshaping after the VPM, in which the 
original pupil is restored. IH] 
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Figure 5: The intensity pattern in the reimaged pupil plane beyond a VPM, 
for an illustrative telescope containing three off-axis mirrors. 

2 The non-ideal PSF — partial AO correction 

A coronagraph on a ground-based telescope will probably be designed to deal 
with a partially-corrected image, where the effects of atmospheric turbulence 
have been partially removed. Such designs have been made, for example, for 
the VLT [ISl [IS] and Gemini [T7] telescopes. 

In what follows I will assume that the central obstruction either does not 
exist or has been removed by Guyon's method. As Figure H] shows, this is 
important. 

The Point Spread Function (PSF) of an adaptive optics system is of 
course affected by a large number of parameters, but when residual phase 
errors are small the PSF is well approximated by a diffraction-limited core, 
flanked by a halo whose functional form is given by the power spectrum of 
the residual phase. An optimized AO system can largely remove low spatial 
frequencies from the atmospheric phase power spectrum, but scales within 
or comparable to actuator spacings cannot be removed. The overall effect is 
that the diffraction-limited part of the PSF sits in a "hole" in the power-law 
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halo. 

Some analytical progress is possible in the case of well-corrected images, 
and this provides a valuable guide to what may be expected and a check in 
various cases of numerical simulations. This is useful in coronagraphic appli- 
cations where a large dynamic range is enquired in the calculation and small 
numerical errors may become relatively important. The essential results are 
the effect of a PVM on a partially-corrected image, and some knowledge of 
the throughput of these systems for targets that are not quite on-axis (the 
targets of interest of course). 

It is possible to calculate the form of the PSF in the limiting case of 
good correction with high actuator density, for a coronagraph built around a 
VPM. The calculation is outlined in Appendix 1, and shows that the coherent 
diffraction-limited core is simply removed from the PSF, with only the halo 
(the residual phase power spectrum) remaining. This is intuitively what one 
would expect, since the VPM is fundamentally an interference device. 

In Appendix 2 it is shown that the throughput of a OVC, for an off- 
axis source, is a strong power-law function of offset distance, with the index 
depending on the charge of the mask. 

Finally, in Appendix 3 it is shown that a wide variety of useful phase 
masks are, in effect, superpositions of VP Ms. 

Taken together, these results suggest that the attainable dynamic range of 
OVCs and related systems, in the partially-corrected AO case, will converge 
at large radii to values fixed by the AO system itself, while more marked 
variations in performance will be apparent at small radii. To investigate this 
in more detail, three general-purpose coronagraphic designs will be studied: 
OVCs for charges 2 and 6, a quadrant phase mask [9], and a prolate-apodized 
mask [TS]. The prolate mask affects amplitude, not phase, and so gives 
excellent performance at the expense of a loss of light. The example used 
here is the A = 0.99 case reported by Soummer et al., which transmits only 
26% of the incident energy. 

Calculating the performance of these systems in the case of partial AO 
correction requires Monte Carlo simulation. The technique used here for 
numerical simulations is discussed in detail by Sivaramakrishnan et a/..|19]. 
These authors describe a useful way of summarizing the AO system perfor- 
mance with just two parameters, the initial uncorrected turbulence strength 
(the usual D/ro) and an effective actuator spacing d^Q. Two illustrative 
AO systems are modelled here to investigate the performance of the OVC. 
The first is a fairly high-performance AO system, with D/r^ = 10 and 
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D/d^Q = 25. Here the diffraction-limited part of the PSF is quite large 
(extending to about 8\/D. The second models a more general-purpose AO 
system, with D/tq = 10 and D/d^Q = 10. 

Some care has to be taken in these numerical simulations as very large 
arrays would be needed both to eliminate low-level aliasing effects and also to 
give adequate detail in the computed PSFs. For the results to be discussed, 
an aperture of 128 pixels diameter was embedded in a 1024 square array. The 
results were checked explicitly against analytical results, both the absolute 
accuracy and for trends. 

Aliasing is an issue for the high levels of accuracy that are required in this 
type of simulation. The VPM diffracts light symmetrically outside the pupil 
according to a power law, and this slow decline aliases back into the pupil 
where there should be zero amplitude. The effect is small but perceptible. 
It can be largely eliminated by deriving the amplitude within the pupil for 
a diffraction-limited source, and then subtracting this from calculations in 
other cases. One check on this method is to compare analytical solutions 
for an obstructed telescope with the simulations corrected in this way (see 
Figure E]) - the method works well for the VPMs, giving reliable results at 
the 10^ level of dynamic range. The quadrant mask simulations work well 
to at least this level without any corrections, probably because of the strong 
angular modulation of extra-pupil diffracted light (Figured]). 

The results are shown in Figures [7| and [HI The diffraction-limited "core" 
region is apparent, transitioning into the power-law "halo" . The halo is unaf- 
fected by the coronagraphs and the interesting region is within the "shoulder" 
that marks the transition. The extent of this region is determined by the 
turbulence strength and the actuator density. 
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Figure 6: A comparison of analytical results for the PSF with numerical 
calculations for a diffraction-limited telescope with a 20% linear central ob- 
struction. Colours code the no- vortex case m = (black), charges 2 (red) 
and 6 (green). Dashed lines are for the numerical calculations, which are for 
a 128-pixel diameter aperture embedded in a 1024-pixel square array, with 
an aliasing correction as described in the text. 
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Figure 7: Simulation results for the long-exposure point spread function of an 
on-axis source, for a the modestly-corrected AO case (no central obstruction, 
(igg = ro). The intensities are normalized to the diffraction-limited case. 
Colours code the no- vortex case m = (black), charges 2 (red) and 6 (green), 
as well as the quadrant phase mask (orange) and prolate-apodized mask 
(blue). Dashed lines show the relative throughput for an off-axis source, as a 
function of off-axis distance. Calculations for the quadrant phase mask, here 
and elsewhere, were at the optimum 45° to the lines of phase discontinuity. 
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Figure 8: Simulation results for the long-exposure point spread function 
of an on-axis source, for a well-corrected AO case (no central obstruction, 
2.5(igg = To). The intensities are normalized to the diffraction-limited case. 
Colours code the no- vortex case m = (black), charges 2 (red) and 6 (green), 
as well as the quadrant phase mask (orange) and prolate-apodized mask 
(blue). Dashed lines show the relative throughput for an off-axis source, as 
a function of off-axis distance. 
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Each coronagraphic method works similarly outside a few \/D. As ex- 
pected, the coherent core of the PSF is largely removed, and so the attainable 
dynamic range is simply set by the strength of the PSF halo at the shoul- 
der. There are strong differences in effectiveness within ~ 2\/D. These are 
ambitiously small offsets to consider, since aberrations and scattering in the 
telescope are extremely difficult to control at these levels. Nonetheless, the 
ultimate purpose of coronagraphs is to be effective at a few \/D. 

The question then arises of how well the various coronagraphs would 
transmit an off-axis source. What is the exclusion radius Tc beyond which a 
source is unaffected by the coronagraph? 

In the case of the OVC can be estimated by requiring that the phase 
introduced by the OVC (charge m), across the off-axis source, should be less 
than about tt. Hence 

m A 

~ 

suggesting that high-charge OVCs will not be able to detect companion 
sources as close as low-charge ones. In fact, it is possible to calculate an- 
alytically the response of an OVC to off-axis sources in the diffraction limit 
(see Appendix 2) and this shows that, close to on-axis, the response is a very 
strong function of charge. Larger charges attenuate very strongly close to 
the axis. Combining these two results suggests that overall higher charges 
should have bigger exclusion radii and is a useful analytical verification of 
the trends in the simulation results. 

For the quadrant mask a narrow radius of exclusion is expected, roughly 
Vc = {2/'k){\/ D) and rather like a charge 2 OVC. The prolate case is too 
complicated to estimate because of diffraction at a stop in the focal plane. 

As a source emerges "from the vortex" the PSF is quite distorted (Figure 
[H]). To quantify the exclusion radius I used as a metric the flux through a 
small circular aperture, extending out to the first (diffraction limited) dark 
ring and centered at the input position of the off-axis test source. This is 
normalized to the flux through the same aperture for the partially-corrected 
PSF, with no coronagraph. The results for the four coronagraphs are super- 
imposed on Figures [3 and [HI 
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Figure 9: A gray-scale image illustrating a source emerging from the vortex 
in the m = 6 well=corrected case. 
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Figure 10: The attainable contrast through an aperture, as defined in the 
text. The analytical results would suggest that performance would converge 
(at least for the phase masks) at large radii and low turbulence, as observed 
in these results. 
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Figure [TUl plots the dynamic range, defined as the flux through the aper- 
ture of an off-axis source, divided by the flux though this aperture from the 
off-axis source, as attenuated by the coronagraph. The general trends are as 
expected; all designs converge to a similar value at large radii, where perfor- 
mance is flxed by AO parameters. Inwards of this, as would be expected, the 
charge 6 design does somewhat better than charge 2 until a point is reached 
where the wider exclusion radius of the charge 6 mask attenuates the off- 
axis source markedly. These trends are rather similar for both levels of AO 
correction, although of course the attainable dynamic range is much better 
for a better-corrected wavefront. What is most apparent however is that 
the quadrant mask does better than the two pure vortex masks, and seems 
to have a well-balanced mixture of various charges to perform well at both 
large and small radii. The performance however will depend on angle relative 
to the phase discontinuities in the mask, which produce "blind lines". The 
prolate-apodized mask also performs well if the loss of light is not an issue. 

A charge-2 OVC would be a good choice for an instrument which had no 
angular dependence of its response in the focal plane, but is more complex to 
manufacture than a quadrant phase mask. Whether there is an "optimum" 
mask of this type in some sense remains an open question. 

3 Conclusions 

The phase vortex mask is an interesting basis for a coronagraph because it 
offers the possibility of complete on-axis extinction of a pure Airy diffraction 
pattern. This result can be demonstrated analytically for even topological 
charges. Moreover, masks of even charges form the Fourier basis to build 
up many other varieties of phase mask, for example the quadrant mask, all 
of which will share the extinction property, and analytical tractability, of a 
single-charge mask. Odd charges always degrade performance so the class of 
useful masks is restricted by this fact. 

A ground-based telescope however will typically have two aspects that 
degrade the performance of vortex phase masks, and phase masks derived 
from them. The central obstruction of the telescope, and the effects of at- 
mospheric turbulence (even partially corrected) are limitations. "Removing" 
the central obstruction, by Guyon's method for example, seems essential to 
attain anything more than modest dynamic range, a conclusion which ap- 
plies to either diffraction-limited images or moderate levels of AO correction. 
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This of course is a common problem for coronagraphs. 

The effect of this class of phase masks can be analyzed quite generally 
when partially corrected turbulence is present, and the conclusion is that 
they will tend to remove the coherent core of the AO-corrected images and 
leave the incoherent power-law halo. The implied outer working radius is 
the "shoulder" in a typical AO image where the core meets the halo, largely 
independently of the details of the mask but fixed by details of the AO system 
such as the actuator density. 

The inner working radius is however affected by the details of the mask. 
It is determined by the extinction of a hypothetical source adjacent to the 
on-axis one. Here the masks do differ, and it is possible to show that high- 
charge masks have a wider zone of extinction than low-charge ones. Taken 
together with the residual light that is transmitted from the on-axis source, 
the detailed simulations show that there is rather little to choose between the 
various masks considered here. In fact the quadrant mask is as good as any. 
High-charge masks seem to offer relatively small gains, for the particular 
turbulence parameters that were examined. The quadrant mask has the 
advantage of being relatively simple to manufacture |2U] . by comparison with 
a vortex where a smooth gradient in phase has to be achieved over the mask. 
Overall however the phase masks are not the limiting factor in dealing with 
AO imagery and it seems likely that a suitable mask could be designed for 
most applications by suitable tailoring of the basis set of vortex masks. 

Appendix 1: Effect of the vortex filter on a 
partially- corrected image 

In the Fourier optics approximation, as is well known, the mathematical 
description of a coronagraph requires three Fourier transforms. The first 
operates in the exit pupil of the telescope or adaptive optics system, and 
produces the amplitude distribution in the focal plane. This amplitude is 
modified by a stop or phase plate. A further (inverse) transform gives the 
resultant amplitude in a reimaged pupil plane, and here a Lyot or pupil stop 
is typically used to eliminate unwanted diffraction. A third transform is used 
to obtain the amplitude at a final focal plane. I will denote co-ordinates in 
the first pupil by r, the first focal plane by k, in the second pupil plane by 
r' and in the final focal plane hj k'. A circular, unobstructed aperture is 
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assumed. 

For a well-corrcctcd adaptive optics system, the amplitude in the exit 
pupil is given by exp(t(f){r)) ~ 1 + i(f){r). If the number of actuators is large 
(the mathematical requirement is that the scale of correlation of the phase 
is much shorter than the pupil size) the a standard argument gives the long- 
exposure point spread function (PSF). This is the diffraction-limited PSF of 
the pupil, scaled by the Strehl ratio, plus a halo which is proportional to 
the power spectrum of the phase, scaled by (1-Strchl). (The appearance of 
a second-order quantity, namely the power spectrum of the phase, indicates 
that for consistency the phase exponential has to be expanded to second 
order. Carry these second-order terms through the calculation gives the 
appropriate scalings of "core" and "halo" by phase variance or Strehl ratio.) 

Subject to the same assumptions, it is possible to find the PSF in the final 
focal plane, after application of a vortex phase filter in the first focal plane, 
and a pupil stop in the reimaged pupil. For an unobstructed diffraction- 
limited system, the final PSF is of course zero; for partial correction, the 
PSF is just the phase power spectrum. The vortex filter removes the coherent 
component, the Airy disc, and only the incoherent part passes through. 

I now give the details of this argument, which, although cumbersome, 
involves only the repeated application of the physical approximations already 
stated. For simplicity, constants of proportionality that are unnecessary to 
the main argument will be omitted. 

In the first focal plane, the amplitude will be (proportional to is under- 
stood from here on) - 

A{k) = A{k) + 1 j dr(f){r) exp{ik ■ r)Q(r) 

in which k denotes a (suitably scaled) position in the focal plane, A is 
the Airy amplitude pattern, and fl{r) is a function which is unity inside the 
pupil and zero outside. 

This amplitude is multiplied by a vortex filter function f{k) which is of 
the form exp(zm^), m being an even integer and 9 being the azimuthal angle 
in the focal plane. //* = 1 will be a useful property later in the argument. 

Propagating to the reimaged focal plane, and applying the full-size pupil 
stop there gives an amplitude leaving that pupil which is 

n{r') [ dkA{k)f{k)exp(-ik-r'). 
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Here r' is the co-ordinate in the reimaged pupiL The extinction property 
of the even-charge vortex filter is that 

n(r^) j dkA{k)f{k) ex.p{-ik ■ r') = 0, 

so the amphtude leaving the reimaged pupil is 

iVL{r') j dk (^j dr (p{r) exp{zk ■ f^Q^f^^ f{k) exp{—zk ■ r'). 

Finally, therefore, the amplitude in the second focal plane will be 

A'{k') = j dr' iQ:{r') j dk j dr (p{f^exp{ik-f^il{f^f{k)exp{—ik-r')exp{ik'-r'). 

The time- or ensemble- averaged intensity is 

I{k') = {A'{k')A'*{k')] 
To carry out the sixfold integration and averaging, write 

A'*{k') = J dr[ (-z)n(r{) J dki J dn </'(rl) exp(-zA^i • n)n{ri) f {ki)* x 
X exp(iki ■ r'l) exp(— z/c^ • r'l). 

in which the subscript 1 is used to indicate corresponding dummy vari- 
ables of integration. 

Re-arranging the order of integration, the first part of the sixfold integra- 
tion involves r and fl : 

J J drdri {(f){f^(f){fi)) exp (ik ■ f — iki ■ r^ij fl{f)Q{fi). 

Writing fl = f+ u, the integral contains the term 

((/)(r)(/)(f + u)) ~ Ciu) 

which is approximately equal to the spatially-invariant phase autocorre- 
lation function C (the equality cannot be exact because of the edges of the 
pupil) . If C declines rapidly with u (many actuators) , then the u and f in- 
tegrals separate to give the approximate result for this first pair of integrals 
as 
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in which P is the phase power spectrum. This will turn out to be all that 

passes through the vortex coronagraph. 
The next pair of integrals is 

dkdkj{k)f{kiy exp (^-ik ■ ? + ik^ ■ rl) V{ki)A{k - h). 

The Airy amplitude A is strongly peaked compared to other terms, and 
picks out k = ki. The phase function / has unit modulus and so disappears, 
leaving 

J dkexp (^—ik ■ r' + ik ■ r[j V{k) — C{r'i — fl). 

The final pair of integrals is 

J J drdri exp (^-ik' ■ ? - ik' ■ r]) C(ri - fi)9.{?)n{P^) . 

With the same change of variable, and the same argument about the scale 
of C compared to the pupil size, it follows that 

I{k') ~ P{k') 

so that the effect of the vortex filter is to remove the coherent Airy core of 
the PSF and leave only the incoherent halo. 




Appendix 2: Analytical estimate of off-axis be- 
haviour of a vortex 

In the high-Strchl approximation the optical vortex filter acts mainly on the 
diffraction-limited core of the AO-corrected PSF. It is important to know 
what happens when this core is not centered on the vortex filter, as this 
will be the case for point sources of interest near to a bright star. In the 
high-Strehl limit, where the amplitude is linear in the residual phase, the 
filter will act on the off-centered Airy disc but pass the uncorrelated (halo) 
light, as is the case for an on-axis source. In this section I show that the 
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transmitted energy is a simple power-law function of the offset distance in 
the focal plane. 

Let us therefore consider an Airy pattern, off-axis by an amount s at some 
angle a (with respect to the origin of the focal plane co-ordinate system {k, 6). 
To simplify the discussion, I work in scaled units of length in both the pupil 
plane and focal plane. In the pupil plane, scale lengths by the pupil radius, 
and in the focal plane, scale by the factor FX/tt, F being the focal ratio of 
the system and A the operating wavelength. This means that the Fourier 
transforms and PSFs can be written compactly. 

The functional form of the off-axis Airy pattern is the standard one in a 
variable k, with 

K = + - 2ks cos(a - 9)) 

k being the scaled radial co-ordinate in the focal plane system centered 
on the vortex filter. The Airy pattern is therefore given (to within constants 
irrelevant to this discussion) by Ji(k)/k. 

To write this in the needed variables, namely k and s, I apply the Gegen- 
bauer identityPTj. Applying this to the present case, 

= 2 2^(l + ri) Q(cos(a-^)) (4) 

K n=0 ^ ^ 

in which C\ is a Gegenbauer polynomial [22]. 

To evaluate the transmission of an off-axis source through the centered 
vortex filter, we can perform a Fourier transform, term by term, of the am- 
plitude after the vortex filter 

exp(«mt^j 

in which the transform is with respect to the variables k and ^, and carries 
the intensity into a re-imaged focal plane in which the co-ordinates are r. 

If the source is on axis, there is no amplitude within the re-imaged pupil 
f2(r) - this is the coronagraphic behaviour of the vortex filter. Off- axis, we 
can find the leading order term in s for which there is a finite amplitude 
within fi. This gives an estimate of the near-to-axis transmission behaviour 
of the vortex filter. 

Off axis (non-zero s) is more complicated mathematically but involves no 
new physics. The steps are as follows. 
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(1) In the Fourier transform into the re- imaged pupil 



J dk exp tmO ^ exp (tk ■ -rj 



we note that 

k ■ r = kr cos(/5 — 6) 

in which the polar co-ordinates in the re-imaged pupil are r and (3. The 
transform is therefore 

J kdk J d6 exp{tm6) ^ exp {tkrcos{f3 — 6)) 

(2) We then insert the Gegenbauer expansion for Ji{n)/n. Expanding 
the Gegenbauer polynomials, and using multiple-angle formulas, rather than 
powers, we find that the angular integral in the transform involves terms 
which are of the form 

dO exp {imO) exp {ikr cos(/3 — 6)) cos(p(a — 0)) 

with m and p both integers. These integrals can be systematically eval- 
uated by using the identity [23] 

exp -(t--))= t'^Uz). (5) 

^ n=-oo 

with the substitution t = exp(z(/3 — 9)), followed by multiplication by 
the appropriate cosine factor and then term-by-term integration. The result 
involves terms in Jp-m and Jp+mi as well as trigonometrical factors in (3. 

(3) The purpose of this calculation is to evaluate the transmission of the 
coronagraph at small s, so at this stage we limit the Gegenbauer series by 
expanding the Jn+i{s)/s term to third order in s. This limits the expansion 
to fourth order Bessel functions and a third order Gegenbauer polynomial. 

(4) We now have the fc-integral to do. This resolves into a number of 
integrals of the Weber-Schafheitlin type (pTj, namely 

J^{at)J,{bt) 
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which have a discontinuity at a = 6 and different behaviour at a > 6 and 
a < b. In our case, this behaviour corresponds to whether any hght is 
diffracted beyond the vortex filter into the re- imaged pupiL 

(5) We carried out this lengthy series of calculations in the computer 
algebra package MATHEMATICA, and investigated topological charges on 
the vortex masks of m = 2, 3, 4, 5, 6. What is of interest is to note the first 
power of s, the focal-plane offset, that gives any diffracted light within the re- 
imaged pupil. If this power is 7 then this means that the transmission energy 
of an off-centered image is proportional to (s/sq)'^'^. For odd topological 
charges, light is diffracted into the pupil regardless of offset; this corresponds 
to the known result that these charges are ineffective in a coronagraph. For 
even charges, we find the interesting result that 27 = m, so that the exponent 
equals the charge. This means that high-charge vortex filters will have a much 
steeper response to off-axis objects. 

(6) Finally, we evaluate sq by integrating the leading-order term in s over 
the re-imaged pupil. What is of interest here is the relative behaviour, so 
we evaluate with respect to the energy passed by a m = filter. This takes 
care of the various multiplicative constants that have been ignored so far. In 
(focal plane) scaled variables, the transmission of energy I for charges 2, 4 
and 6 is given by 



h 
h 
h 



6 

32 
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This result is illustrated in Figure [TT] 
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0.6 0.7 



offset / FX/n 

Figure 11: The transmitted intensity through a vortex, as a function of 
off-axis angle, for charges 2 (red), 4 (blue) and 6 (green). The analytical 
approximation only applies at small s and to emphasize this, lines are shown 
dashed beyond s = 0.5 
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Here is a specific example of tlie calculation for a charge 2 vortex. Ex- 
panding the Gegenbauer series (Equation H]) to third order in s gives 



K k 



^ ^ f Uk) cos{e-a) \ ^ 



2 , -3Ji(A;) + 3J3(A;) + QJ^{k) cos(20 - 2a)\ 



—2J2{k) cos{6 — a) + Ji{k) cos{6 — a) + J^^ik) cos(36' — Sa) 

2Ak 



The Fourier transform back into the pupil involves the term-by-term in- 
tegration 

J kdkjde exp{z29) exp {ikrcos{l3 - 9)) 

in which we specialize to the charge 2 case as an illustration. The variable 
r is the scaled co-ordinate in the re-imaged pupil plane. 

The first term (no dependence on offset s) for the pupil-plane amplitude 
A is 

A{r) = < r < 1 

-27rexp(2«/5) 



r > 1. 



The next term (linear in s) is 

A{r) = znsexp^i^a + P)) r 0<r<l 

—insexp{—i{a — 3/3)) 



r > 1. 



Here we see the first appearance of finite amplitude within the pupil. 
For larger charges the calculations have to be carried to higher in s before 
this happens. Figure [12] illustrates the appearance of the reimaged pupil, 
calculated in this analytical fashion, for a small offset of a source in the case 
of a charge 6 vortex. 
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Figure 12: The transmitted intensity (stretched scale) in the pupil through a 
charge 6 vortex, for a small offset of the source in the focal plane. This was 
calculated analytically by the method in the text. 

A useful check is to repeat this calculation for no vortex, namely a charge 
of zero. Although complicated expressions result, they simplify as expected; 
the intensity is uniform over the pupil, with a constant phase gradient. 

Appendix 3. The quadrant mask as a super- 
position of vortex masks 

The quadrant mask amplitude is defined as a function of azimuthal angle 9 
in the focal plane as 
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AO) = 


1 0<e <n/2 


AO) = 


-1 7r/2<e<7T 


Ai9) = 


1 n <9 < 3n/2 


AiO) = 


-1 3n/2 <e <27i 



with no radial dependence. This amphtude can be expressed as a Fourier 
series: 

A{e) = -(... + -e-i°^^ + -e~''' + e-'^' + e^'' + -e^^' + -e^o*^ + . . .) (6) 

so, a sum of various evenly charged vortex phase masks. 

The effect of this (or other phase masks represented as Fourier series) can 
be analyzed analytically for an Airy diffraction pattern. Diffraction from the 
focal plane will involve integrals of the form 

/ drr-^ de A{9) exp{ikr cosie - p)) 
Jo k Jo 

which requires the result 

/ d9 exp{in9) exp{tz cos{9 — /3)) = 2tt exp{ml3)Jn{z) 
Jo 

for z real and n an integer. (This can be deduced from Equation [5] by the 
method described there.) 

The radial integrals are then of the Weber-Schafheitlin form, so the am- 
plitude in the reimaged pupil is a weighted sum of the amplitudes for the 
various vortex masks that appear in Equation [61 Each of these is zero within 
the pupil as long as only evenly charged vortices appear in Equation [61 

As an example. Figure [T3l shows a numerical and an analytical calculation 
of the reimaged pupil plane after application of a quadrant phase mask. This 
technique allows analytical calculations to be performed for any useful phase 
mask which does not have a radial dependence of the phase. 
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Figure 13: The transmitted intensity (stretched scale) in the pupil through 
a quadrant phase mask, calculated numerically (right) and analytically with 
ten non-zero terms in the Fourier expansion. 
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